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Evidence for a Bound H-dibaryon from Lattice QCD 
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We present evidence for the existence of a bound H-dibaryon, an / = 0, J = 0, s = —2 state with 
valence quark structure uuddss, at a pion mass of ~ 389 MeV. Using the results of Lattice 
QCD calculations performed on four ensembles of anisotropic clover gauge-field configurations, with 
spatial extents of L ~ 2.0, 2.5, 3.0 and 3.9 fm at a spatial lattice spacing of bs ~ 0.123 fm, we find 
an H-dibaryon bound by B^ = 16.6 ± 2.1 ± 4.6 MeV at a pion mass of ~ 389 MeV. 



It is now well established that quantum chromodynam- 
ics (QCD), the theory describing the dynamics of quarks 
and gluons, and the electroweak interactions, underlie 
all of nuclear physics, from the hadronic mass spec- 
trum to the synthesis of heavy elements in stars. To 
date, there have been few quantitative connections be- 
tween nuclear physics and QCD, but fortunately. Lat- 
tice QCD is entering an era in which precise predictions 
for hadronic quantities with quantifiable errors are being 
made. This development is particularly important for 
processes which are difficult to explore in the laboratory, 
such as hyperon-hyperon and hyperon-nucleon interac- 
tions for which knowledge is scarce, primarily due to the 
short lifetimes of the hyperons, but which may impact the 
late-stages of supernovae evolution. In this letter we re- 
port strong evidence for a bound H-dibaryon, a six-quark 
hadron with valence structure uuddss, from n/ = 2 -f 1 
Lattice QCD calculations at light-quark masses that give 
the pion a mass of ttItt ^ 389 MeV. 

The prediction of a relatively deeply bound system 
with the quantum numbers of AA (called the H-dibaryon) 
by Jaffe [1] in the late 1970s, based upon a bag-model 
calculation, started a vigorous search for such a sys- 
tem, both experimentally and also with alternate the- 
oretical tools. Experimental constraints on, and phe- 
nomenological models of, the H-dibaryon can be found 
in Refs. [SlISlll]. While experimental studies of doubly- 
strange hypernuclei restrict the H-dibaryon to be un- 
bound or to have a small binding energy, the most recent 



constraints on the existence of the H-dibaryon come from 
heavy- ion collisions at RHIC, from which it is concluded 
that the H-dibaryon does not exist in the mass region 
2.136 < Mh < 2.231 GeV [5], effectively ehminating the 
possibility of a loosely-bound H-dibaryon at the physical 
light-quark masses. Recent experiments at KEK suggest 
there is a resonance near threshold in the H-dibaryon 
channel (S]. 

The first study of baryon-baryon interactions with Lat- 
tice QCD was performed more than a decade ago [Zl [5]. 
This calculation was quenched and with rriTr^ 550 MeV. 
The NPLQCD collaboration performed the first fif — 2 + 
1 QCD calculations of baryon-baryon interactions [9l [10] 
at low-energies but at unphysical pion masses. Quenched 
and dynamical calculations were subsequently performed 
by the HALQCD collaboration [III dl]. A number of 
quenched Lattice QCD calculations [H [H [HJ [H HH HH] 
have searched for the H-dibaryon, but to date no defini- 
tive results have been reported. Earlier work concluded 
that the H-dibaryon does not exist as a stable hadron 
in quenched QCD [T7] , while more recent work [THl US] 
finds a hint of a bound state. By inserting energy- and 
sink-dependent potentials into the Schrodinger equation 
in the SU(3) limit, a hint of an H-dibaryon has been 
found in Ref. |20j . however, this hint evaporates when 
SU(3)-breaking is included [Hj. 

In this work, Liischer's method [22l |23l [24l [25] is em- 
ployed to extract two-particle scattering amplitudes be- 
low inelastic thresholds from Lattice QCD calculations. 
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In the situation where only a single scattering channel is 
kinematically allowed, the deviation of the energy eigen- 
values of the two-hadron system in the lattice volume 
from the sum of the single-hadron masses is related to the 
scattering phase shift, 6{q), as is made explicit in eq. ([T]). 
The Euclidean time behavior of Lattice QCD correlation 
functions of the form C^{t) = (0|x(t)x'f (0)|0), where x 
represents an interpolating operator with the quantum 
numbers of the one-particle or two-particle systems un- 
der consideration, determines the ground state energies 
of the one-particle and two-particle systems, Ei = m 
and E2, respectively (we focus only on the ground state 
of the two-particle system in this work). The form of 
the interpolating operators, and the methodology used 
for extracting the energy shift are discussed in detail in 
Ref. [55]. For gauge- field configurations that have dif- 
ferent lattice spacings in the temporal and spatial di- 
rections (anisotropic lattices), the two-particle energy is 
given by E2 = 2-^ g^/C^ + rn?, where ^ — hg/ht is the 
lattice anisotropy. By computing the mass of the parti- 
cle and the ground-state energy of the two-particle sys- 
tem, the squared momentum, (in spatial lattice units 
(s.l.u)), which can be either positive or negative, is de- 
termined by this relation. For s-wave scattering below 
inelastic thresholds, is related to the real part of the in- 
verse scattering amplitude through the eigenvalue equa- 
tion [23] (neglecting phase-shifts in Z > 4 partial-waves) 

.cot%) = -L.(,^(A)^^, (1) 

where the S-function is given by 

|j|<A , 

S{x) = lim y 47r A . (2) 

^ ' A^oo ^ |j|2 -X ^ ' 

This relation provides a Lattice QCD determination of 
the value of the phase shift at the momentum \/q^. 

Determining energy-levels with the same quantum 
numbers in multiple volumes allows for bound states to 
be distinguished from scattering states. A bound state 
corresponds to a pole in the S matrix, and in the case of 
a single scattering channel, is signaled by cot (5(g) +i 
in the large volume limit. Writing q = in for two-particle 
states that are negatively shifted in energy, E2 < 2m, in 
the lattice volume, the volume dependence of the bind- 
ing momentum in the large volume limit follows directly 
from eq. ([!]) and is of the form 25] 

n = 7 + |( e-^^ + V2 e-^^^ )+..., (3) 

where 7 is the infinite-volume value of the binding mo- 
mentum, under the assumption that 7 <C m^, and 171 
is treated as a fit parameter. With calculations in two 
or more lattice volumes that both have < and 
qcot5{q) < it is possible to perform an extrapolation 



with eq. ([3]) to the infinite-volume limit to determine the 
binding energy of the bound state. Boo = 7^/to. The 
range of nuclear interactions is set by the pion mass, 
and therefore the use of Liischer's method requires that 
niT^L ^ 1 in order to strongly suppress the contributions 
that depend upon the volume as e^"^"^ p7] . 

Our present results are from calculations on four en- 
sembles of n/ = 2-1-1 anisotropic clover gauge- field config- 
urations at a pion mass of ^ 389 MeV, a spatial lat- 
tice spacing of bg ^ 0.123(1) fm, an anisotropy [221153] of 
^ = 3.50(3), and with spatial-extents of 16, 20, 24, 32 lat- 
tice sites, corresponding to spatial dimensions of L ~ 2.0, 
2.5, 3.0 and 3.9 fm respectively, and temporal extents of 
128, 128, 128, and 256 lattice sites, respectively. The pre- 
cision of the calculations is sufficiently high that the ex- 
ponential volume dependence of the single baryon masses 
can be cleanly quantified. The A mass, unlike that of 
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FIG. 1: Left panel: the mass of the A as a function of e"™" 
where L is the spatial extent of the lattice. The left-most 
(red) point and uncertainty is the infinite-volume extrapola- 
tion of the other (blue) points calculated in lattice volumes 
with spatial extents of, from left-to-right, L = 32,24,20, and 
16. The curve corresponds to the best straight-line fit. Right 
panel: the energy-momentum relation of the A calculated on 
the 32^ X 256 ensemble. The points (and uncertainties) are 
the results of lattice calculations and the (red) curve corre- 
sponds to the best fit (see text). The units of the vertical axes 
in both plots are t.l.u., and of the horizontal axis of the right 
plot are (t.l.u.)^ 

the TT and kaon, is found to have statistically signifi- 
cant volume-dependence, as shown in the left panel of 
fig. [1] It is clear that the A mass on the 16"^ x 128 en- 
semble (mTfi ~ 3.9) is significantly higher than its 
infinite-volume value and, more importantly, is shifted 
by an amount that is comparable to the two-baryon en- 
ergy shifts. The deviation found in calculations on the 
20^ X 128 ensemble (m^ri — 4.8) is much less than that 
of the 16^ X 128 ensemble, but we choose to use only cal- 
culations on the 24"^ x 128 ensemble {rriT^L = 5.8) and on 
the 32^ X 256 ensemble [mT^L ~ 7.7) in the bound-state 
analysis. 

Liischer's method assumes that the continuum single- 
hadron energy-momentum relation is satisfied over the 
range of energies used in the eigenvalue equation in 
eq. ([!]). In order to verify that this is the case, sin- 
gle hadron correlation functions were formed with well- 
defined lattice spatial momentum, k = for |np < 5. 
As the low-lying states in the lattice volume have energies 
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that are small compared with the A mass, it is sufficient 
to determine the non-relativistic energy-momentum rela- 
tion, 



2Mi 



8M| 



(4) 



The A energy as a function of momentum calculated 
on the 32'^ x 256 ensemble is shown in the right panel 
of fig. [ij and yields Mo,Mi,M2 of 0.22135(10)(05), 
0.2231(34)(13), 0.261(26)(04) t.l.u, respectively. Clearly 
the special-relativity limit of Mq = Mi = M2 is satisfied, 
but an uncertainty of ^ 2% is introduced into from 
the uncertainties in the energy-momentum relation. The 
use of relativistic or lattice dispersion relations leads to 
similar conclusions. 

In the absence of interactions, the AA-S7V-ES system 
is expected to exhibit three low-lying states as the mass- 
splittings between the single-particle states are (on the 
32^ X 256 ensemble) 



2(Ms-Ma) = 0.01317(13)(19) t.l.u 
+ Mn-2Ma ^ 0.003397(61)(65) t.l.u . (5) 



However, if interactions generate a bound state, it is ex- 
pected that the splitting between the ground-state and 
the two additional states will be larger than estimates 
based upon the single baryon rest masses. The effective 
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FIG. 2: The EMPs for the A (left panels) and the splitting 
between the AA system and twice the A mass (right panels) 
calculated on the 24^ x 128 (upper) and 32^ x 256 (lower) 
ensembles. The units of both axes are t.l.u. 

mass plot (EMP) for the A calculated on the 24^ x 128 
and 32'^ x 256 ensembles that have been optimized for the 
ground-states using the matrix-Prony method jl9j are 
shown in the left panels of fig. [2] and clear plateaus are 
identified. The calculated EMP for the energy-splittings 
between the AA-SA^-SS coupled-channels (optimized us- 
ing the matrix-Prony method) and twice the energy of 
the A (formed from the ratio of correlation functions) on 
the 24'^ X 128 and 32^^ x 256 ensembles are shown in the 
right panels of fig- [2] The finite- volume binding ener- 
gies calculated in the L = 16, 20, 24 and 32 lattice vol- 
umes are 12.3(1.1)(4.0), 4.5(1. 1)(1.3), 16.3(1. 2)(1.4), and 



16.6(1.4)(3.1) MeV, respectively. In each lattice volume, 
the results are consistent with a single isolated ground- 
state with an energy that is below the AA threshold (and 
considerably below the EN and EE thresholds). The 
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FIG. 3: The results of the Lattice QCD calculations of 
—icotS versus q'^/ml- obtained using eq. ([l]), along with the 
infinite- volume extrapolation using eq. ([3|. The dark (blue) 
(light (green)) lines correspond to the statistical (systematic 
and statistical uncertainties combined in quadrature) 68% 
confidence intervals calculated on the 24"^ x 128 ensemble 
(lower) and 32'^ x 256 (upper) ensembles. The (red) point and 
its uncertainty at —icotS = -1-1 corresponds to the infinite- 
volume extrapolation, the inner uncertainty being statistical 
and the outer being the systematic and statistical combined 
in quadrature. 

energy splittings and their uncertainties extracted from 
both ensembles lead to negative values of q cot S indicat- 
ing that they both lie on the bound-state branch of the 
S-function (eq. ([2|), and thus leads us to identify the 
H-dibaryon. The extracted values of —i cot 6 from the 
24^ X 128 and 32^ x 256 ensembles and their uncertain- 
ties are shown in fig. [3] along with the infinite-volume 
extrapolation implicit in eq. ([T]), and made explicit in 
eq. ^ . The H-dibaryon binding energy at this pion mass 
is found to be 



16.6 ± 2.1 ± 4.5 ± 1.0 ± 0.6 MeV 



(6) 



where the first uncertainty is statistical, the second sys- 
tematic, the third is an estimate of the uncertainty in 
the infinite-volume extrapolation, and the fourth is the 
uncertainty from the energy-momentum relation. Com- 
bining the various systematic uncertainties in quadrature 
gives = 16.6 ±2.1 ±4.6 MeV. A Monte-Carlo propa- 
gation of the uncertainties indicates that there is a prob- 
ability greater than 0.98 that the H-dibaryon is bound at 
this pion mass. 

In conclusion, we have presented strong evidence for 
the existence of a bound H-dibaryon from Lattice QCD 
calculations at a pion mass of ~ 389 MeV. Our cal- 
culations were performed in four lattice volumes, and a 
negatively shifted ground-state was found in all four vol- 
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umes. In order to avoid contamination from finite- volume 
modifications to the A mass and interactions, only the re- 
sults obtained in the larger two volumes were used in the 
infinite-volume extrapolation. Within the uncertainties, 
the ground-state energies in the largest two volumes are 
the same, indicating that both volumes are large com- 
pared with the H-dibaryon size. This is consistent with 
the calculated binding energy. Calculations were per- 
formed at only one lattice spacing. However, given that 
lattice-spacing artifacts in these calculations are expected 
to scale as 0(6^), we expect such contributions to be 
small. Moreover, general arguments based on the low- 
energy effective theory of the Symanzik action suggest 
that 0{b^) effects largely cancel in forming the energy 
difference. Consequently, we expect the observation of 
the H-dibaryon to survive the continuum extrapolation. 
However, the quark-mass dependence of the H-dibaryon 
binding energy is presently unknown, so a direct com- 
parison of our result with experiment is not yet possible. 
As with all such lattice calculations, we cannot rule out 
the possibility of an additional deeper bound state of 
the same quantum numbers in this channel that couples 
weakly to the interpolating operators. 

The results of the Lattice QCD calculations presented 
in this letter provide the first clear evidence for a bound- 
state of two baryons directly from QCD. This is further 
strong motivation for pursuing Lattice QCD calculations 
in larger volumes, at smaller lattice spacings, and over 
a range of light-quark masses including those of nature, 
as the present calculations demonstrate that the study of 
light (hyper-) nuclei directly from QCD is feasible. 
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